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Abstract 

It has been claimed (Canup and Ward 2002; Ward 2003) that a long-lived 
massive (compared to the mass of the Galilean satellites) circumplanetary gas 
disk is inconsistent with Jupiter's low obliquity. Such a constraint could be 
downplayed on the basis that it deals with a single observation. Here we argue 
that this argument is flawed because it assumes a solar system much like that 
of the present day with the one exception of a circumjovian disk which is then 
allowed to dissipate on a long timescale (10 6 — 10 7 yrs). Given that the sequence 
of events in solar-system history that fit known constraints is non-unique, we 
choose for the sake of clarity of exposition the orbital architecture framework 
of Tsiganis et al. (2005), in which Jupiter and Saturn were once in closer, less 
inclined orbits than they are at present, and show that Jupiter's low obliquity is 
consistent with the SEMM (solids-enhanced minimum mass) satellite formation 
model of Mosqueira and Estrada (2003a,b). 

1 Introduction 

While there may be good reason to believe that planet-disk interactions may excite the 
eccentricities of at least isolated planets (Goldreich and Sari 2003; Ogilvie and Lubow 
2003) there are currently no studies to guide our understanding of likely outcomes in 
the case of multiple planet (or satellite) systems. Therefore, one may consider circular, 
coplanar giant planets as a starting condition. Tsiganis et al. (2005) have recently 
argued that evolution through the 1:2 Jupiter-Saturn mean motion resonance (MMR) 
of a quasi-circular, coplanar and compact solar system that is allowed to evolve by 
planetesimal scattering is consistent with the observed semi-major axes, eccentricities 
and mutual inclinations of Jupiter, Saturn, Uranus and Neptune. In a companion 
publication, Gomes et al. (2005) argue that this resonance crossing may be linked 
to the Late Heavy Bombardment of the terrestrial planets taking place ~ 700 Myr 
after their formation (Hartmann et al. 2000). For such a scenario to work the bulk of 
the divergent migration between Jupiter and Saturn and the resonance passage itself 
must take place following gas dissipation, and also later than simulations with an 
evenly spread disk of planetesimals would indicate. Whether or not the identification 
with the Late Heavy Bombardment is correct, in this viewpoint the solar system 
must have been more compact and regular prior to gas dissipation (during its first 
10 6 — 10 7 years) than it is today. The authors also briefly consider the potentially 
disruptive consequences of such a scenario for the regular and irregular satellites of the 



giant planets, but conclude that at least the regular satellites might have "survived" 
unscathed. 1 

Jupiter's low obliquity ~ 3° may be indicative of its formation by hydrodynamic 
gas accretion. Yet, it has been noted that secular spin-orbit resonances can compli- 
cate this straightforward interpretation. In particular, adiabatic passage through a 
resonance matching the spin axis precession rate to the v\§ precession frequency of the 
orbit plane due to the gravitational perturbation of Saturn (Hamilton and Ward 2002; 
Canup and Ward 2002; Ward 2003) and u 17 due to the gravitational perturbation of 
Uranus (Hamilton, pers. comm.) may result in obliquities significantly larger than 
observed. For Jupiter the amplitude of the z/ 16 term (with a period of ~ 50, 000 years) 
is ~ 0°.36 and the i>n term (with a period of P±q ~ 450,000 yrs) is ~ 0°.055, which 
could have resulted in obliquities of up to ~ 26° and ~ 14°, respectively, had these 
resonances been crossed adiabatically as the circumplanetary gas disk was dissipated 
(Ward 2003). This leads Canup and Ward (2002) to argue that the circumplanetary 
disk must have viscously evolved in a timescale sufficiently short (O(10 5 ) yr) as to 
preclude adiabatic passage, resulting in a gas-starved (or gas-poor [Mosqueira et al. 
2000; Estrada and Mosqueira 2005]) satellite disk. 

The issue we tackle here is whether our decaying turbulence 2 SEMM model 
(Mosqueira and Estrada 2003a,b; hereafter MEa,b) is especially susceptible to secu- 
lar spin-orbit resonances, and inconsistent with Jupiter's low obliquity. In particular, 
we focus on the v\q term. The reasons for this are: First, the period of the orbital 
precession vyj [Pyj ~ 450, 000 yrs) is already very close to the precession period of 
Jupiter's spin axis due to the solar torque on the Galilean satellites (which are locked 
to the Jupiter's equator plane by this planet's oblateness; Goldreich 1965). A slight 
adjustment of satellite or planet positions might be enough to place Jupiter spin axis 
precession period in one side or the other of this resonance, so that any formation 
model is apt to be affected. Second, even if the resonance is crossed the limiting 
obliquity for adiabatic passage is significantly smaller than that for the v\§ term. 
Furthermore, the timescale for adiabatic passage in the case of the vyj secular spin- 
orbit resonance may be longer (O(10 7 ) yrs) than the timescale for gas dissipation by 
photoevaporation 3 . Third, it is likely that Uranus and Neptune formed after Jupiter 
and Saturn, once most of gas in the planetary and circumplanetary disks had already 
dissipated. Thus, from here on we focus on v\q and consider Jupiter and Saturn only. 

1 Though the issue is not spelled out in detail (presumably due to space limitations), in the case 
of the Saturnian satellite system Iapetus' relatively low eccentricity e ~ 0.03 may be the main 
constraint for such a scenario. For this reason, it would appear unlikely that Titan could owe its 
eccentricity to a close encounter between Saturn and another giant planet. 

2 Consistent with numerical simulations that show turbulence decay in the absence of a "stirring" 
mechanism (Hawley et al. 1999). 

3 Shu et al. (1993) conclude that EUV from the central star may photoevaporate a T Tauri disk 
in ~ 10 7 yrs outside of a gravitational radius r g ~ 10 AU. But there is considerable uncertainty in 
this estimate. For instance, Adams et al. (2004) argue that the disk is likely to be photoevaporated 
by FUV radiation from neighboring stars. Furthermore, these authors state that photoevaporation 
is effective to a significantly smaller radius 0.1 — 0.2r g . 
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Figure 1: Inclinations of Jupiter and Saturn in degrees derived from a two-body 
secular perturbation theory calculated over a time span of 2 x 10 5 yrs. (a) Case 
centered on 1983 using actual orbital elements (Murray and Dermott 1999). (b) Case 
in which Jupiter has semimajor axis aj = 5.45 AU, Saturn has as = 8.5 AU, and 
their mutual inclination is Ijs = 10 -3 radians. 

2 Secular Perturbation Theory 

In a solar system consisting of Jupiter and Saturn it is straightforward to construct 
a Laplace-Lagrange secular solution (e.g., Brouwer and Clemence 1961; Murray and 
Dermott 1999). For / << 1, the orbital inclination I and the longitude of the 
ascending node with respect to the invariant plane Q are given by 

I sin Q = 7 1 sin7 1 + I 2 sm(f 2 t + 72), (1) 

and 

IcosQ = /1COS71 + I 2 cos(f 2 t + 72), (2) 

where f 2 is the eigenfrequency, l\ and I 2 are eigenvector components, and 71 and 
72 are phases. If we use parameters for Jupiter and Saturn as observed, we obtain 
f 2 = —7.06 x 10~ 3 ° yr -1 and I 2 = —6.30 x 10~ 3 (in radians). This yields a secular 
oscillation with period of Piq ~ 51,000 yr. These values are not too dissimilar from 
the secular solution of the planetary system (e.g., Murray and Dermott 1999). 

Let us now consider a time early on before Jupiter and Saturn had passed through 
the 1:2 MMR but after most of the planetary gas disk had dissipated by photoevap- 
oration in a timescale 10 6 — 10 7 yr. Since the planetary nebula may shield it to some 




degree, it may be appropriate to assume that the much denser subnebula takes longer 
to dissipate. If so, at this time both the precession of Jupiter's spin axis and orbital 
plane may be significantly faster than they are today. Because of the scattering on 
nearby planetesimals (Gomes et al. 2005) and the dissipation of the nebula, Saturn 
would dominate the precession of Jupiter's orbital plane. 4 We take as nominal values 
for the orbital parameters of Jupiter and Saturn the starting conditions in Tsiganis 
et al. (2005), namely, aj = 5.45 AU, as = 8.50 AU (a few tenths of an AU inside 
the 1:2 MMR) and sin I js ~ 10~ 3 , where Ijs is the relative inclination of the orbital 
planes of Jupiter and Saturn. With these parameters 5 we obtain f 2 = —1.60 x 10~ 2 ° 
yr _1 and I 2 = —2.72 x 10~" 4 (in radians). As shown in Fig. 1, in this case the period 
of Jupiter's orbital plane precession is P±e ~ 23, 000 yrs. 



3 Obliquity Variation 

Jupiter's current spin axis precession period is ~ 4.5 x 10 5 yrs due mostly to the 
solar torque exerted on the Galilean satellites (e.g., Ward 1975) 6 . However, a massive 
circumplanetary gas disk would result in a much shorter precession period. The 
inner parts of the disk (r < r t = (2M p /M Q J 2 R 2 p af) 1/b AORj, where M p = Mj, 
a p = 5.45 AU, R p ~ 2Rj, and J 2 oc 1/R P ~ 0.008 are the planetary mass, semi-major 
axis, radius and quadrupole gravitational harmonic following envelope collapse [which 
may take place in a fast 10 4 — 10 5 yrs timescale, Hubickyj O., pers. comm.], and M Q 
and Rj are the Sun's mass and Jupiter's present radius) would orbit in the plane of 
the planet's equator and precess as a unit with the planet, whereas the outer parts of 
the disk (r > r t ) would not (Goldreich 1966). The spin axis s then precesses around 
the orbit normal n at a rate given by (e.g., Tremaine 1991) 



— = a(s-n)(sx n), (3) 

where s and n are unit vectors, and the precession constant 7 is given by 

4 Ward (1981) treats the two-orbit/nebula problem as dispersal takes place. Here we consider the 
case when the nebula (but not the subnebula) has already been dissipated by some means. 

5 An estimate for the gap opening timescale may be obtained using the tidal torque formula (Lin 
and Papaloizou 1993) T gap ~ (A/a) 5 P/fi 2 , where A is the gap's half-width, fj, is the mass ratio of the 
secondary to the primary, and P is the orbital period of the secondary. Using A ~ (a s ~aj)/2 ~ 1.5 
AU, we find that Jupiter and Saturn would have cleared the gas disk in between in a short timescale 
t ~ 10 4 yrs. In particular, this timescale is shorter than satellite formation timescale 10 4 — 10 6 
(depending on location; Mosqueira et al. 2001) in MEa,b. 

6 We obtain a precession period of 4.8 x 10 5 yrs using a moment of inertia for Jupiter of K = 0.26, 
which yields a spin angular momentum J p = 4.4 x 10 45 g cm 2 /s. There is some uncertainty in the 
moment of inertia ~ 5% (Fortney J., pers. comm.) but this doesn't affect the argument. What is 
important to note, however, is that about 33% of the precession constant a is contributed by the 
torque of the Sun directly on the planet because of its oblateness. Adding a circumplanetary gas 
disk decreases this fraction, which we then ignore. 

7 In actuality, there should be another term added to this precession rate due to the torque of 



<x=~2jf I Z(a)a 3 da, (4) 
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where the surface density of circumplanetary disk E(a) is assumed to drop-off sharply 
at r D ~ 2r c , where r c = Rh/A& ~ 15Rj is the centrifugal radius (MEa), fl p is the 
planet's orbital frequency and the total angular momentum of the precessing system 
is given by 



rD 



H = J P + 2n(GM p ) l/2 / E(a)a 3/2 da, (5) 

J Rp 

where G is the gravitational constant and J p is the spin angular momentum of the 
planet. For E oc 1/a, we can write a = M D Vt 2 r 2 D / (AH) and H = J p +2/3M D ^GM p r D , 
where M D is the disk mass. The precession period is given by T = 2n/(a cos 9), where 
cos 9 — s • n is the obliquity. Given that in our SEMM model M D ~ 10M sats , where 
M sa ts ~ 4 x 10 26 g is the mass of the Galilean satellites 8 , the spin of the planet 
provides ~ 90% of the angular momentum of the system H ~ 5 x 10 45 g cm 2 /s (see 
MEa Table 3). Using an obliquity of 9 ~ 3° and t\d = AORj (which implies a surface 
density E~2x 10 4 g/cm 2 at 15Rj consistent with the value obtained by applying 
the inviscid gap-opening criterion to Ganymede in a disk with aspect ratio ~ 0.1 
[MEb]), we obtain a precession period of T ~ 4 x 10 4 yrs, which is slightly shorter 
than Jupiter's current orbital plane secular precession period P 16 ~5x 10 4 yrs, but 
it is longer than Jupiter's orbital plane precession period when Saturn was inside the 
1:2 MMR, i.e., P w ~ 2 x 10 4 yrs. Hence, in our decaying turbulence SEMM model 
it is possible for Jupiter either to have crossed the secular spin-orbit resonance before 
the Keplerian disk reached its quiescent phase (and satellites formed) at a time when 
the viscous evolution of the disk was likely driven by Roche-lobe gas inflow and the 
resonance passage was non-adiabatic, or not to have crossed this resonance at all. 
However, it may be possible to alter this conclusion by choosing different parameters, 
such as a more massive and extended circumplanetary disk. Thus, next we consider 
the case of resonance passage. 

The limiting obliquity 9 max that could be generated by the obliquity "kick" in- 
curred during adiabatic resonance passage in the non-capture direction (a < 0) (which 
implies resonance crossing as the circumplanetary gas disk is dissipated and the spin- 
axis precession period increases) is given by (Henrard and Murigande 1987; Ward and 
Hamilton 2004) 

the extended part of the disk r > r t on the inner part of disk. However, the contribution of this 
extended region drops rapidly with distance, and in our model the gas surface density drops-off at 
a radial location ro < rt. 

8 There is some ambiguity here because in the model of MEa,b a significant fraction of the mass 
of Callisto is derived from the extended part of the disk, whereas Mo is the mass of the inner disk 
out to about Callisto. On the other hand, Io and Europa should be reconstituted for unaccreted 
volatiles. 
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which yields 9 

max ~ 26° using the current value for I/2I ~ 0°.36, but a significantly 
smaller value of 9 max ~ 9.1° using \I 2 \ ~ 0°.015 obtained from our nominal, low 
mutual inclination case. Furthermore, the minimum time for adiabatic crossing is 
(Ward and Hamilton 2004) 



Taking 9 ~ 9°, we find ~ 2 x 10 8 yrs for our nominal case, which is much longer 
than the gas dissipation timescale. This means that at least for the nominal case 
the crossing must be non-adiabatic and the final obliquity is rate dependent. Taking 
9 ~ 3°, we calculate r min ~ 9 x 10 4 yrs for current solar-system parameters, but it 
is T min ~ 2 x 10 7 yrs (here P 16 ~ 23, 000 yrs is about a factor of two shorter than 
the present value, but this is more than compensated by the much smaller value of 
1^2 1) f° r our nominal case, which is comparable to or longer than the gas dissipation 
timescale. We may estimate the resulting obliquity to be 9 < tan 1 / 3 (|/ 2 |) < 4° for 
our nominal case, which is consistent with Jupiter's observed value. At any rate, this 
would all be irrelevant if the resonance were not crossed. 

4 Conclusions 

We have briefly investigated the consequences for Jupiter's obliquity of a satellite 
formation model in which a massive (compared to the Galilean satellites yet enhanced 
in solids by a factor of ~ 10 compared to the solar composition minimum mass model) 
subnebula is allowed to dissipate in a long timescale (10 6 — 10 7 yrs) after the dispersal 
of the nebula itself. For the sake of specificity, we have adopted the model of Tsiganis 
et al. (2005) in which the solar system was more compact and regular before Jupiter 
and Saturn crossed the 1:2 MMR than it is today. We find that such a combined 
scenario (by no means unique) does not imply the likelihood of a larger obliquity for 
Jupiter than is observed. This is both because the secular z/ 16 spin-orbit resonance may 
not be crossed, and because the resulting obliquity may be consistent with Jupiter's 
value even if it is crossed. We conclude that Jupiter's low obliquity is compatible with 
our SEMM satellite formation model (MEa,b) provided one allows for solar-system 
conditions early-on unlike those presently observed. 
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